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We propose a method to generate analytical quantum Bell inequalities based on the principle 
of Macroscopic Locality. By imposing locality over binary processings of virtual macroscopic in- 
tensities, we establish a correspondence between Bell inequalities and quantum Bell inequalities in 
bipartite scenarios with dichotomic observables. We discuss how to improve the latter approximation 
and how to extend our ideas to scenarios with more than two outcomes per setting. 
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I. INTRODUCTION 

Any probability distribution achievable in the mea- 
surement of a bipartite quantum system can be written 

as 



P(a,b\x,y) = Tr(pU x a 



where p is a density operator of unrestricted dimension 
and the ITs are projectors. The set of quantum cor- 
relations (or quantum set) is the set of all probability 
distributions that can be written in this form. It has 
long been known that some correlations in this set are 
nonlocal, in the sense that they cannot be produced by 
pre-established agreement [1]. On another front, the 
quantum set is strictly smaller than the no-signaling 
set, i.e. those correlations that cannot be used to send 
a signal Q. 

The mathematical characterization of the quantum 
set is not trivial [3| . One of the few definite results was 
obtained for the case of two dichotomic measurements 
per site, i.e., x, y, a and b can take only two values 
each (a situation that we shall refer to as 2222). In this 
scenario, up to symmetries, the only Bell inequality is 
the Clauser-Horne-Shimony-Holt (CHSH) Q 



J2(-l) xy (a x b y ) < 2, 



(1) 



where for convenience of notation we have chosen 
x,y G {0,1}, a, 6 G { — 1,+1}. Tsirelson and Lan- 
dau (and, much later, Masanes) independently proved 
that the set of quantum correlations is bounded by a 
quantum Bell inequality obtained from CHSH upon re- 
placing (a x b y ) by ^ a,Ycs'm(a x b y ) [El-Q- More recently, 
Navascues, Pironio and Acm proved that, if we relax 
the tensor structure assumption in (pQ) and instead de- 
mand commutativity between Alice's and Bob's mea- 
surement operators, the resulting set of correlations 



can be approximated, and ultimately reached, by a hi- 
erarchy of semi-definite programs [1, In the pro- 
cess, they proved that the replacement (a x b y ) — > 



^ arcsin D xy , where 



D 



(a x b y ) - (a x }(a y ) 



xy 



(by) 2 ) 



(2) 



in the CHSH inequality leads to a quantum Bell in- 
equality that is stronger than the previously known 
one and identifies the set denoted Q 1 , the first step of 
the semi-definite hierarchy. 

Aside from from the mathematical line of attack, 
several authors have followed the spirit of the Popescu- 
Rohrlich paper [2[ and tried to identify physical prin- 
ciples that would define the quantum set fl0l-[l4{ . The 
principle of information causality has already proved 
useful to derive bounds on the quantum set in non- 
trivial scenarios [Hj]. In this paper, we use the princi- 
ple of macroscopic locality [13] to derive quantum Bell 
inequalities. The recipe is applicable to all bipartite 
scenarios. Focusing on nm22 scenarios, we prove that 
a quantum Bell inequality can be obtained from any 
standard Bell inequality through (j2j). Such inequalities 
are the only analytical approximations to the quantum 
set available to date. 



II. TOOLS 
A. Macroscopic locality 

Before presenting our results, we review briefly the 
idea of macroscopic locality (ML) [l3|. We consider a 
bipartite system, with a source emitting pairs of parti- 
cles to Alice and Bob. Alice can choose a measurement 
x G X = {1 . . . , tua} from a set of tua possible set- 
tings, each producing oIa possible outcomes, denoted 
as a G A = {1,...,gU}- Similarly, Bob can choose a 
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measurement y from the set Y = {l,...,ra#}, with 
ds outcomes b G 5 = {1, . . . , de}. One calls a micro- 
scopic experiment one in which single-particle detec- 
tions are possible. In this case, routine sampling leads 
to reconstructing P(a,b\x,y). By contrast, a macro- 
scopic experiment is defined as follows: in each run, 
N pairs are sent; Alice measures all her particles with 
setting x and records the number n a of particles that 
produced the outcome a G A; Bob does similarly. Af- 
ter repeating this procedure several times, one can re- 
construct P{n a, ft B \x, y) , where n A = [n a =i, n a=dA ] 
and Hb = [nb=i, —<,nb=d B ]- Clearly, some microscopic 
information is lost as soon as N > 1: Alice and Bob do 
not know which of their outcomes came from the same 
pair. 

The principle of ML demands that, in the limit N — >• 
oo and under the assumption that Alice's and Bob's 
measurement devices cannot detect particle number 
fluctuations smaller than O(V^V), the macroscopic 
statistics should not violate any Bell inequality. The 
principle is not trivial: for instance, one can easily 
check that, if the source were a Popescu-Rohrlich (PR) 
box, the macroscopic statistics would violate CHSH up 
to 4, just as the microscopic statistic do. In fact, ML 
is respected if and only if the microscopic statistics 
P(a, b\x, y) belong to Q 1 



B. Recipe for quantum Bell inequalities 

It is now easy to formulate our recipe to generate 
quantum Bell inequalities. Consider a data processing 
of the macroscopic variables 



n A 



a G A', 



(3) 



This defines a scenario with txia measurements for Al- 
ice, vfiB for Bob, \A!\ = d! A outcomes for Alice and 
\B'\ = d' B for Bob - in short, denoted as mArriBd f A d B 
scenario. From ML it follows that, if P(a,b\x,y) be- 
longs to the quantum set, P(a, /3\x,y) cannot violate 
any Bell inequality. In other words, for any Bell in- 
equality 1 < in the scenario mAmBd f A d B , the con- 
dition l[P(a, f3\x, y)\ < defines a quantum Bell in- 
equality for the microscopic correlations. If in addition 
P(a,fi\x,y) can be expressed as a function of the pa- 
rameters of P(a,b\x,y), an analytical bound on the 
quantum set is obtained. 

Of course, we know already that none of these crite- 
ria will come closer to the quantum set than Q 1 does; 
and in general, the bound may even be less tight than 
Q 1 (we show an explicit example below). However, to 
date, this is the only known systematic approach to ob- 
tain analytical approximations to the quantum set in 
any scenario. In what follows, we study some specific 



examples of quantum Bell inequalities generated with 
this recipe. 



III. CASE STUDIES 
A. Two outcomes: Sign-binning 

Let dA = ds = 2 and for convenience we use 
A = B = { — 1,1}. Consider first a specific choice of 
measurements x, y (so, for simplicity, we omit them in 
the notation). From her macroscopic statistics, Alice 
can infer the microscopic marginal average 

(a) = (An A )/N with An A = n a = + i - n a= _i . (4) 

A natural data processing that gives back binary out- 
comes {A! = { — 1,1}) consists in comparing An a with 
(a) in each run: 



Sign-binning: a 



{ +1 if 
\ -1 ol 



AnA > (a) 



otherwise 



(5) 



Bob can of course do the same. Now we want to re- 
late a and f3 obtained by sign-binning to the properties 
of the microscopic distribution. Any macroscopic run 
uses N pairs: denote therefore {o' 1 ',...^'^} the re- 
sults Alice would have observed on those TV pairs in 
a microscopic experiment, and similarly {b^\ b^} 
for Bob. Then we define 



N 



a = sgn 



- (a) \ _ 



V/c=l 



N 



sgn(a / ), (6) 



^ = s g n (e &( % w )^ sgn(&/) - (7) 

In the limit N oo, due to the central limit theorem, 
the joint probability distribution of the variables (a 7 , b') 
is a gaussian distribution Gr(a f , b') with zero mean and 
covariance matrix, Y determined by the microscopic 
distribution [13| 



I -(a) 2 (ab)-(a)(b) 
(ab)-(a)(b) 1 -(b) 2 



(8) 



Then, with the notation / = 
Jda'db' Gp(a / , b')f(a', b f ), a straightforward inte- 
gration shows that 



a = f3 = 0, a/3 



■ arcsinD^ . 



(9) 



All that we have done is valid for any pair of mea- 
surements x, y. Therefore, as claimed, the replacement 
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([2]) is not specific to the CHSH scenario: given any 
bipartite Bell inequality with binary outputs, one can 
impose that P(a, fi\x,y) respects it and then replace 
all the a^r, /3 X and a x f3 y with the corresponding expres- 
sions above: the resulting inequality will automatically 
define a quantum Bell inequality for the microscopic 
distribution P(a, b\x,y). Let us denote the set of corre- 
lations satisfying these inequalities as Q SB . Of course, 
from the above discussion it follows that Q 1 C Q SB . 

Many bipartite Bell inequalities with binary out- 
put are known [H, [l7| , some of which have intriguing 
properties (e.g., they are maximally violated by non- 
maximally entangled states, or define dimension wit- 
nesses). For the sake of this paper, we consider only 
two examples in what follows. 



B. Sign binning for 2n22 

With CHSH being a special case (n = 2) of this 
scenario, we note that the only relevant Bell inequ ality 
in 2n22, up to symmetries, is the CHSH inequality jl8| . 
Using the replacement in (J9j), we obtain the analytical 
inequality 

| arcsin P^ + arcsinZ^i + arcsin Pij — arcsinD2j| < tt 

(10) 

for alH, j = 1, . . . , n as a quantum Bell inequality. For 
the 2222 scenario, we recover the inequality given in 
pjl, thus proving that Q SB = Q 1 in that case. In fact, 
Q 1 = Q SB also holds in any 2n22 scenario: this is not 
unexpected, since it is known that CHSH is the only 
relevant inequality for those scenarios [18j. We give 
an explicit proof in the Appendix [A] In summary, in 
2n22 scenarios, the distribution P{nA,nB\%,y) is local 
if and only if P(a, /3\x,y) is local. This remarkable 
coincidence is however not general as shown in the next 
paragraph. 



C. 3322 Scenario: sign-binning and beyond 

For 3322, besides CHSH, there is only one inequiva- 
lent Bell inequality denoted as 13322 HH G3 • The set 
Q SB is therefore defined by both (fTDj) and 



arcsin D xy — arcsin P32 — arcsin D 

x+y<4: 



23 



< 2tt. 
(11) 



By considering a specific slice of the 3322 no-signaling 
poly tope, we found that this set of inequalities gives 
weaker constraints than ML, i.e. Q 1 C Q SB '. This can 



be seen from Figure [TJ In the two dimensional slice 
defined by the three extremal points, Pi, P2 and P n , 
qSB giveg a j3 0un d of 7 3 322 = 0.4 while ML (i.e. Q 1 ) 
gives 0.2. 
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FIG. 1: A slice of the 3322 no-signaling polytope denned by 
the two extremal points Pi, P2 (13322 = 1) and P n (random 
noise). In this slice, the bound for the local set L and Q are 
the same, which is 0. The notation of the points is taken 
from M~SI] and the extremality of Pi and P2 follows from 



It is not astonishing that the quantum Bell inequal- 
ities obtained from sign binning may be weaker. Geo- 
metrically, sign binning is obtained by slicing the Gaus- 
sian curve into two symmetrical parts, each part rep- 
resenting a macroscopic outcome. In general though, 
macroscopic locality is recovered by slicing the curve 
into infinitely many parts or macroscopic outcomes 
pjl , as shown in Figure [3l 

This motivates us to investigate different forms of 
binning. Here we show a form of data processing in- 
duced by the Gaussian function 



P(x)= exp(-^). 
The binning is then defined as follows, 



(12) 



with probability P(Sn) 
3-binning: a = ^ +1 otherwise and Sn > 
- 1 otherwise 



(13) 



where Sn = An a — (a). This generates a new set of 
probabilities {P(c/, /3'\x, y)} in the 3333-scenario, all of 
which can be computed analytically. Thus, this allows 
us, in principle, to have an analytical quantum Bell 
inequality. Here, we numerically determine the bound 
and find that the new set, Q SB gives a better approxi- 
mation to Q. This can be seen in Figure [T] where Q SB 
gives 0.304 by taking a = 0.028 in Eq. ([12]). 
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D. Binning for n > 2 microscopic oul :omes. 

The sign binning processing is of course a natural 
processing for binary outcome settings, though not op- 
timal. It can also be applied to settings with more 
than two outcomes; but here we sketch a more natu- 
ral processing for such scenarios. As an example, we 
will focus on the case of 2233: X = Y = {0, 1} and 
A = B = {0,1,2}. We first define the quantity 



rii(a) -pfN 



N 



,Vi eA 



(14) 



for a particular fixed value of x and y. The data pro- 
cessing, in this case called the triangle-binning, is then 
defined as 



a 



/c , s . t . Iu 



(15) 



Bob can similarly do the same. In the limit N — >• oo, 
the probability distribution for the variables (a x ,/3 y ) 
is a Gaussian distribution Gr(& x ,Py) with the covari- 
ance matrix T defined in the usual sense. The for- 
mulas for the probabilities P(a, /3\x,y) include definite 
trigonometric integrals; they may not lead to elegant 
analytical expressions, but they can then be evalu- 
ated numerically to check the method. It is known 
that 2233 has only one relevant Bell inequality of the 
Collins-Gisin-Linden-Massar-Popescu (CGLMP) type 
[HI: 1 cglmp < 0. We can then impose the condition 
Icglmp[P(o£, f3\x, y)\ < 0. The result is shown in Fig- 
ure [2] where the constrained set of correlation is now 
denoted as Q TB . Contrary to sign binning, triangle 
binning does not reproduce the predictions of ML. Tri- 
angle binning may thus not be worth pursuing further, 
but it shows that generalization is indeed possible. 



n 



L Q Q Ppr 

-r-H 1 



0.5 0.73 0.75 1 

FIG. 2: The line of random marginals denned by the ran- 
dom noise, P n and the generalized PR box, Ppr [23j]. The 
bounds given by the local set, L, macroscopic locality, Q 1 



and triangle binning, 
tively. 



Q are 0.5, 0.75 and 0.73 respec- 



IV. CONCLUSION 

In this paper, we have studied a new post-processing 
method, called sign binning, to define a set of quantum 
correlations denoted as Q SB , which necessarily con- 
tains the quantum set, Q. Sign binning, as mentioned, 




FIG. 3: a) ML in the sense of binning into infinitely many 
parts, b) Sign binning for scenarios with binary outcomes, 
c) 3-binning induced by the probability, P(x). d) Triangle- 
binning for scenarios with three outcomes. 



cannot give a bound on the quantum set that is tighter 
than Q 1 , but it does give a straightforward, intuitive 
way of generating quantum Bell inequalities. To our 
knowledge, this recipe is the only method to date to 
derive analytical quantum Bell inequalities. 
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Appendix A: Locality of Sign Binning is 
Equivalent to Macroscopic Locality in 2n22 

Lemma 1. Let T be an n + 2 square matrix of the form 



r = 



A C 
C T B 



(Al) 



where C is a given 2 x n real matrix and A, B are such 
that An = Bjj = 1 for i = 1, 2 and j = 1, n. Then 
there exists a choice of the remaining entries of A and 
B such that T > iff there exists x G [—1,1] such that 

l-x 2 - Cl - Cl + 2xC u C 2l > 0, (A2) 

fori = l,...,n. 

Proof Let us first prove that, if condition (|A2|) holds, 
then T can be made positive semidefinite. Suppose 
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that, indeed, such an x exists and 
can take A to be 



x\ < 1. Then, we 



A = 



1 x 
x 1 



> 0. 



(A3) 



According to Schur's theorem [24|, if A > 0, a matrix 
of the form (|Aip is positive semidefinite iff B' = B — 
C T A~ X C > 0. Since the non-diagonal entries of B are 
not determined a priori, we can always choose them 
such that B[- = for i ^ j. To see that B' is positive 
semidefinite, we then only have to show that B f u > 0. 
But 



1 



r<2 r 1 ^ 



2xC u C 2i 



1 



(A4) 



that is non- negative by hypothesis. We have just 
proven that, for \x\ < 1, condition (|A2[) grants posi- 
tive semidefiniteness. Suppose now that (|A2|) holds for 
x = 1. Then the equation reads 



c 2l f > 0, 



for i = 1, n. 



(A5) 



It follows that Cii = for all i. In order to show 
that T can be completed to a positive semidefinite ma- 
trix, take the orthonormal basis {|0), |1)} and define 
the vectors: 



v lj2 = |0); v i+2 = C u \0) + V 1 " ( A6 ) 

Then, the Gram matrix Y[- = Vi-Vj is positive semidef- 
inite, has Is in the diagonal and its off-diagonal sub- 
matrix coincides with C. 

The case x = — 1 can be treated analogously (simply 
take v\ = — #2). 

Now we will prove the opposite implication: suppose 
that there is some way to complete Y such that T > 0. 
Let f be such completion and take x = A\ 2 . If x = ±1, 
then the Gram decomposition of = [24] is such 
that || = ||iT2 1| = 1 and v\ - v 2 — ±1. This implies 
that V\ = d=^2, and so Cu = #1 • v 2 +i = ±^2 ■ = 
d=C2i, and condition ()A2|) holds for x = ±1. 

Suppose that, on the contrary, \x\ < 1. Then A > 0, 
so, by Schur's theorem, > 0, and condition ()A2|) 
holds. 

□ 

Theorem 2. Let Y be a matrix such as the one ap- 
pearing in the definition of the previous lemma. Then, 
T can be made positive semidefinite iff, for all i,j = 
l,...,n ; i ^ j, 

arcsin(Cii) + arcsin(C2i) + arcsin(Cij) — 



arcsin(C2j) 
plus permutations of the minus sign. 



< 7T, (A7) 



Proof. By lemma [TJ positive semidefiniteness is equiv- 
alent to the existence of an x E [— 1, 1] satisfying the 
conditions (|A2|) . Without loss of generality, we assume 
that Cu = sm(<fii),C 2 i = sin(^), for — tt/2 < 0i : <fii < 
7r/2. Then, conditions (|A2|) can be reexpressed as 



- cos 



+ 0i) < x < cos((j)i - Si). 



(A8) 



An x satisfying all these conditions exists iff the min- 
imum of the upper limits is greater than or equal to the 
maximum of the lower limits. In other words, Y can be 
completed iff 



cos 



+ 0j) < cos( 



i),Vi,j. (A9) 



Call a* ee \ pj ee 1^+6^1. Then, < c H ^ 3 < 

7r, and the positivity condition reads 



cos(c^) 



COS 



7) > 0. 



(A10) 



Running through all possibilities ([c^ < 7r/2,/3j < 
tt/2], [a % < tt/2, ft > tt/2], [a, > tt/2, ft < tt/2], [a, > 
7r/2,/3j > 7r/2]), one can check that this condition is 
equivalent to 0.1 + fij < 7r, and so we arrive at equations 
(1X71) . 

□ 

We will now prove the claimed result in the article 
that Q SB = Q 1 for 2n22. This proof is an extension of 
the proof presented in [9]. For the case of 2n22, let Y 1 
be a certificate of order 1 for a particular P(a, b\x,y). 
Then we have 



cl z? 




(All) 



where Xu = Yjj = 1 for i = 1,2, j — l,...,n. Also, 
Ca = (Cai, CU2) and C# = (Cbi, • • • , C^n) are the 
marginal correlations for Alice's and Bob's measure- 
ment respectively. By Schur's Lemma [24], T 1 > is 
equivalent to the positive semidefiniteness of 



X z 
Z T Y 



(A12) 



where the matrix Y has 



{1 - cl„i 



r 2 



diagonal elements 

may assume that all the diagonal elements are non 
zero for n > 2; if any of them is zero, then the 
outcome of that particular measurement in deter- 
ministic and can be accounted for with a local 
hidden variable model. Therefore we can multiply 
T on both sides with the diagonal matrix, M = 



6 



The condition T > is then equivalent to 

r = ( C T B ) " °' (A13) 
where the 2 x n matrix C has elements Cij = (Cij — 



CiC^/^l-Cf^Jl-C] and An = = 1 for % = 

1, 2 and j = 1, n. Now by Lemma [T] and Theorem[2j 
a certificate of order 1 exists, and therefore probability 
distribution, P(a,b\x,y) is Q 1 if and only if condition 
in ([TO]) is satisfied. 
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